Fermi Large Area Telescope Observations of the Cosmic-Ray Induced 7-ray Emission 

of the Earth's Atmosphere 



A. A. Abdo, 1 - 2 M. Ackcrmann, 3 'E M. Ajello, 3 W. B. Atwood, 4 L. Baldini, 5 J. Ballet, 6 G. Barbiellini, 7 - 8 

D. Bastieri, 9 ' 10 B. M. Baughman, 11 K. Bechtol, 3 R. Bellazzini, 5 B. Berenji, 3 E. D. Bloom, 3 E. Bonamente, 12 ' 13 
A. W. Borgland, 3 A. Bouvier, 3 J. Bregeon, 5 A. Brez, 5 M. Brigida, 14 ' 15 P. Bruel, 16 R. Buehler, 3 T. H. Burnett, 17 

S. Buson, 10 G. A. Caliandro, 18 R. A. Cameron, 3 P. A. Caraveo, 19 J. M. Casandjian, 6 C. Cecchi, 12 - 13 
6. Qelik, 20 ' 21 ' 22 E. Charles, 3 A. Chekhtman, 1 ' 23 J. Chiang, 3 S. Ciprini, 13 R. Claus, 3 J. Cohen-Tanugi, 24 
J. Conrad, 25 - 26 ' 27 F. de Palma, 14 ' 15 S. W. Digel, 3 E. do Couto e Silva, 3 P. S. Drell, 3 R. Dubois, 3 D. Dumora, 28 ' 29 
C. Farmer, 24 C. Favuzzi, 14 ' 15 S. J. Fegan, 16 W. B. Focke, 3 P. Fortin, 16 M. Frailis, 30 Y. Fukazawa, 31 S. Funk, 3 -Q 
P. Fusco, 14 ' 15 F. Gargano, 15 N. Gehrels, 20 ' 32 ' 33 S. Germani, 12 ' 13 B. Giebels, 16 N. Giglietto, 14 ' 15 F. Giordano, 14 ' 15 
T. Glanzman, 3 G. Godfrey, 3 I. A. Grenier, 6 M.-H. Grondin, 28 ' 29 J. E. Grove, 1 L. Guillemot, 34 S. Guiriec, 35 

E. Hays, 20 D. Horan, 16 R. E. Hughes, 11 G. Johannesson, 3 A. S. Johnson, 3 T. J. Johnson, 20 ' 33 W. N. Johnson, 1 
T. Kamae, 3 H. Katagiri, 31 J. Kataoka, 36 N. Kawai, 37 ' 38 M. Kerr, 17 J. Knodlseder, 39 M. Kuss, 5 J. Lande, 3 

L. Latronico, 5 M. Lemoine-Goumard, 28 ' 29 F. Longo, 7 ' 8 F. Loparco, 14 ' 15 B. Lott, 28 ' 29 M. N. Lovellette, 1 
P. Lubrano, 12 ' 13 A. Makeev, 1 ' 23 M. N. Mazziotta, 15 J. E. McEnery, 20 ' 33 C. Meurer, 25 ' 26 P. F. Michelson, 3 
W. Mitthumsiri, 3 '[lT. Mizuno, 31 C. Monte, 14 ' 15 M. E. Monzani, 3 A. Morselli, 40 I. V. Moskalenko, 3 S. Murgia, 3 
P. L. Nolan, 3 J. P. Norris, 41 E. Nuss, 24 T. Ohsugi, 31 A. Okumura, 42 N. Omodei, 5 E. Orlando, 43 J. F. Ormes, 41 
D. Paneque, 3 J. H. Panetta, 3 D. Parent, 28 ' 29 V. Pelassa, 24 M. Pepe, 12 ' 13 M. Pesce-Rollins, 5 F. Piron, 24 
T. A. Porter, 4 S. Raino, 14 - 15 R. Rando, 9 - 10 M. Razzano, 5 A. Reimer, 44 ' 3 O. Reimer, 44 - 3 T. Reposeur, 28 ' 29 
L. S. Rochester, 3 A. Y. Rodriguez, 18 M. Roth, 17 H. F.-W. Sadrozinski, 4 A. Sander, 11 P. M. Saz Parkinson 4 
C. Sgro, 5 G. H. Share, 1 ' 45 E. J. Siskind, 46 D. A. Smith, 28 ' 29 P. D. Smith, 11 G. Spandre, 5 P. Spinelli, 14 - 15 
M. S. Strickman, 1 D. J. Suson, 47 H. Takahashi, 31 T. Tanaka, 3 J. B. Thayer, 3 J. G. Thayer, 3 D. J. Thompson, 20 
L. Tibaldo, 9 - 10 ' 6 D. F. Torres, 48 - 18 G. Tosti, 12 ' 13 A. Tramacere, 3 ' 49 Y. Uchiyama, 3 T. L. Usher, 3 V. Vasileiou, 21 ' 22 
N. Vilchez, 39 V. Vitale, 40 - 50 A. P. Waite, 3 P. Wang, 3 B. L. Winer, 11 K. S. Wood, 1 T. Ylinen, 51 ' 52 ' 26 and M. Ziegler 4 

(Fermi-LAT collaboration) 

1 Space Science Division, Naval Research Laboratory, Washington, DC 20375, USA 
2 National Research Council Research Associate, National Academy of Sciences, Washington, DC 20001, USA 

W. W. Hansen Experimental Physics Laboratory, 
Kavli Institute for Particle Astrophysics and Cosmology, 
Department of Physics and SLAC National Accelerator Laboratory, Stanford University, Stanford, CA 94305, USA 
* Santa Cruz Institute for Particle Physics, Department of Physics and Department of Astronomy and Astrophysics, 
University of California at Santa Cruz, Santa Cruz, CA 95064, USA 
5 Istituto Nazionale di Fisica Nucleare, Sezione di Pisa, 1-56127 Pisa, Italy 
6 Laboratoire AIM, CEA-IRFU /CNRS/Universite Paris Diderot, 
Service dAstrophysique, CEA Saclay, 91191 Gif sur Yvette, France 
1 Istituto Nazionale di Fisica Nucleare, Sezione di Trieste, 1-34127 Trieste, Italy 

8 Dipartimento di Fisica, Universita di Trieste, 1-34127 Trieste, Italy 
9 Istituto Nazionale di Fisica Nucleare, Sezione di Padova, 1-35131 Padova, Italy 
10 Dipartimento di Fisica "G. Galilei", Universita di Padova, 1-35131 Padova, Italy 
11 Department of Physics, Center for Cosmology and Astro-Particle Physics, 
The Ohio State University, Columbus, OH 43210, USA 
12 Istituto Nazionale di Fisica Nucleare, Sezione di Perugia, 1-06123 Perugia, Italy 
13 Dipartimento di Fisica, Universita degli Studi di Perugia, 1-06123 Perugia, Italy 
■* Dipartimento di Fisica "M. Merlin" dell' Universita e del Politecnico di Pari, 1-70126 Bari, Italy 
15 Istituto Nazionale di Fisica Nucleare, Sezione di Bari, 70126 Bari, Italy 
16 Laboratoire Leprince-Ringuet, Ecole polytechnique, CNRS/IN2P3, Palaiseau, France 

17 Department of Physics, University of Washington, Seattle, WA 98195-1560, USA 
18 Institut de Ciencies de VEspai (IEEC-CSIC), Campus UAB, 08193 Barcelona, Spain 
19 INAF-Istituto di Astrofisica Spaziale e Fisica Cosmica, 1-20133 Milano, Italy 
20 NASA Goddard Space Flight Center, Greenbelt, MD 20771, USA 
21 Center for Research and Exploration in Space Science and Technology 
(CRESST) and NASA Goddard Space Flight Center, Greenbelt, MD 20771, USA 
22 Department of Physics and Center for Space Sciences and Technology, 
University of Maryland Baltimore County, Baltimore, MD 21250, USA 
23 George Mason University, Fairfax, VA 22030, USA 
2J> Laboratoire de Physique Theorique et Astroparticules, 
Universite Montpellier 2, CNRS/IN2P3, Montpellier, France 



2 



25 Department of Physics, Stockholm University, AlbaNova, SE-106 91 Stockholm, Sweden 
26 The Oskar Klein Centre for Cosmo-particle Physics, AlbaNova, SE-106 91 Stockholm, Sweden 
27 Royal Swedish Academy of Sciences Research Fellow, 
funded by a grant from the K. A. Wallenberg Foundation 
28 Universite de Bordeaux, Centre d' Etudes Nucleaires Bordeaux Gradignan, UMR 5797, Gradignan, 33175, France 
29 CNRS/IN2P3, Centre d' Etudes Nucleaires Bordeaux Gradignan, UMR 5797, Gradignan, 33175, France 
30 Dipartimento di Fisica, Universitd di Udine and Istituto Nazionale di Fisica Nucleare, 
Sezione di Trieste, Gruppo Collegato di Udine, F33100 Udine, Italy 
31 Department of Physical Sciences, Hiroshima University, Higashi-Hiroshima, Hiroshima 739-8526, Japan 
Department of Astronomy and Astrophysics, Pennsylvania State University, University Park, PA 16802, USA 
33 Department of Physics and Department of Astronomy, 
University of Maryland, College Park, MD 20742, USA 
" 14 Max- Planck- Institut fur Radioastronomie, Auf dem Hiigel 69, 53121 Bonn, Germany 
35 Center for Space Plasma and Aeronomic Research (CSPAR), 
University of Alabama in Huntsville, Huntsville, AL 35899, USA 
36 Waseda University, 1-104 Totsukamachi, Shinjuku-ku, Tokyo, 169-8050, Japan 
Department of Physics, Tokyo Institute of Technology, Meguro City, Tokyo 152-8551, Japan 
38 Cosmic Radiation Laboratory, Institute of Physical and Chemical Research (RIKEN), Wako, Saitama 351-0198, Japan 
39 Centre d'Etude Spatiale des Rayonnements, CNRS/UPS, BP 44346, F-30128 Toulouse Cedex 4, France 
4 °Istituto Nazionale di Fisica Nucleare, Sezione di Roma "Tor Vergata", 1-00133 Roma, Italy 
41 Department of Physics and Astronomy, University of Denver, Denver, CO 80208, USA 
4 Department of Physics, Graduate School of Science, 
University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo 113-0033, Japan 
43 Max-Planck Institut fiir extraterrestrische Physik, 85748 Garching, Germany 
44 Institut fiir Astro- und Teilchenphysik and Institut fiir Theoretische Physik, 
Leopold-Franzens-Universitat Innsbruck, A-6020 Innsbruck, Austria 
45 Praxis Inc., Alexandria, VA 22303, USA 
4a NYCB Real-Time Computing Inc., Lattmgtown, NY 11560-1025, USA 
47 Department of Chemistry and Physics, Purdue University Calumet, Hammond, IN 46323-2094, USA 
48 Institucid Catalana de Recerca i Estudis Avancats (ICREA), Barcelona, Spain 
49 Consorzio Interuniversitario per la Fisica Spaziale (CIFS), 1-10133 Torino, Italy 
5 Dipartimento di Fisica, Universitd di Roma "Tor Vergata", 1-00133 Roma, Italy 
51 Department of Physics, Royal Institute of Technology (KTH), AlbaNova, SE-106 91 Stockholm, Sweden 
52 School of Pure and Applied Natural Sciences, University of Kalmar, SE-391 82 Kalmar, Sweden 

(Dated: December 9, 2009) 

We report on measurements of the cosmic-ray induced 7-ray emission of Earth's atmosphere by 
the Large Area Telescope onboard the Fermi Gamma-ray Space Telescope. The LAT has observed 
the Earth during its commissioning phase and with a dedicated Earth-limb following observation 
in September 2008. These measurements yielded ~ 6.4 x 10 6 photons with energies > 100 MeV 
and ~ 250 hours total livetime for the highest quality data selection. This allows the study of the 
spatial and spectral distributions of these photons with unprecedented detail. The spectrum of the 
emission - often referred to as Earth albedo gamma-ray emission - has a power-law shape up to 500 
GeV with spectral index F = 2.79 ± 0.06. 



I. INTRODUCTION 



For an earth-orbiting 7-ray detector such as the Large Area Telescope (LAT) on board the Fermi 7-ray Space 
Telescope (Fermi), the Earth is the brightest source in the sky due to its proximity. The 7-ray emission from the 
Earth is produced by cosmic-ray interactions with the Earth's atmosphere, and is often called the 7-ray albedo. The 
vast majority of cosmic rays are protons and heavier nuclei. Atmospheric 7-ray emission is the result of cosmic- 
ray cascades, mainly through the decay of neutral pions and kaons and through Bremsstrahlung of electrons and 
positrons (dominating at energies below ~50 MeV [l(). Due to the kinematics of the collisions, the cross sections of 
these processes at high energies are peaked in the forward direction. 
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The spectrum of the albedo 7-rays is not uniform across the Earth's disk. The spectrum of the inner part of the disk 
is soft owing to the small number of secondary particles, typically neutral pions, that are back-scattered at large angles 
relative to the direction of the cosmic-ray cascade developing deep into the atmosphere. Cosmic rays that enter the 
atmosphere near grazing incidence produce showers whose forward-moving 7-rays can penetrate the thin atmospheric 
layer making the limb bright when viewed from orbit. The geometric effect of shower production and absorption in 
the atmosphere producing a bright 7-ray horizon was first calculated by Stecker [2L A Monte Carlo model of 7-ray 
production by cosmic-ray interactions in the atmosphere later developed by Morris [l| produced reasonable agreement 
with the spatial and spectral measurements by various balloon-borne and spacecraft instruments. 

The Earth's 7-ray emission has been observed by several previous satellite missions starting as early as the 1960s 
with a detector flown on the third Orbiting Solar Observatory, OSO-3, for energies above 50 MeV This instrument 
detected an enhancement of 7-rays when pointing towards the Earth, consistent with balloon flight measurements 0, 
S B 0) II] and followed up by other satellite missions 0, [l(| HH, 03 • A more detailed study of the 7-ray emission from 
interactions of cosmic rays with the Earth's atmosphere was done by Thompson et al. [13| using approximately 6700 
7-ray events from the Earth, recorded with the second Small Astronomy Satellite, SAS-2. These early measurements 
showed a clear peak of 7-ray intensity towards the Earth horizon, with a factor of ~10 larger intensity than seen 
towards the nadir direction. Also they show a modulation of the 7-ray intensity with respect to the east-west 
direction as expected from the deflection of cosmic rays in the Earth's magnetic field. These findings were confirmed 
and extended to higher energies with the high-energy 7-ray satellite detector on-board the Compton Gamma-Ray 
Observatory (CGRO), the Energetic 7-ray Experiment Telescope (EGRET). EGRET was in operation from 1991 to 
2000, and approximately 60% of the mission's total 5.2 million events were measured from the Earth direction [HI ]. 

The next generation 7-ray observatory, Fermi, was launched on 11 June 2008. The primary instrument on Fermi, 
the Large Area Telescope (LAT), is a pair-conversion telescope that detects 7-rays through conversion into an e + -e~- 
pair. The trajectories of this pair are recorded in the tracker and allows for the reconstruction of the direction of the 
incident photon. The subsequent particle shower development is sampled in the tracker and a calorimeter to determine 
the photon energy. The LAT's wide field of view (~ 2.4 sr), effective area (~ 8000 cm 2 at 1 GeV), improved point 
spread function (PSF), and broad energy coverage (20 MeV to > 300 GeV) provide significantly improved sensitivity 
over its predecessors [15j . 

During routine science operations (the so-called "survey" mode), the LAT points away from the Earth because of 
the relative brightness of the Earth at GeV 7-ray energies and the background it would introduce to astrophysical 
sources. The center of the field-of-view (z-axis) of the LAT is kept within 35° relative to the zenith (the direction 
pointing away from the center of the Earth) for most of the time . However, during the commissioning phase of the 
instrument, the pointing was such that parts of the Earth were within the field-of-view at varying angles with respect 
to the LAT z-axis. In addition, during September 2008 a dedicated observation of the Earth's limb was made. This 
combined dataset can be used for the analysis of the albedo emission. 

There are several reasons that render a study of the Earth's 7-ray albedo with the LAT useful. First, measurements 
of the albedo 7-rays contribute to the understanding of the celestial 7-ray background for the LAT and future 
satellite instruments. Second, such observations can be used for a calibration of the LAT instrument. Finally, if the 
production of secondary 7-rays in the Earth's atmosphere is sufficiently well understood and modeled, observations of 
the albedo emission yield a measurement of the interactions of primary cosmic rays with the Earth's atmosphere and 
can, therefore, provide information on the spectrum of primary cosmic rays, interaction cross sections, atmospheric 
properties, and the deflection of primary cosmic rays by the geomagnetic field. 

The current work describes the analysis of the Earth's 7-rays albedo emission at energies above 100 MeV observed 
with the LAT. 



II. DATASET 

Two main data sets have been used: (i) the Launch & Early Operations (LEO) data taken during the first 60 
days of the mission, during which time the Earth limb was often closer to the field of view than during a regular 
sky survey (only data were used that were suitable for science observations), (ii) two orbits (i.e. close to 3 hours) of 
direct limb observations during which the Earth limb was kept within the field of view with the aim of adding on-axis 
photons to the overall data set. The LEO data contain many different observation modes, but the main modes are a 
north-orbital-pole pointing and a pointing mode with the LAT z-axis centered at the Vela pulsar. All results shown 
here are averaged over many orbits of Fermi and correspond to the highest quality event selection currently available 
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for the LAT data (the so-called "diffuse" class events [lj|). The total number of events > 100 MeV in our dataset is 
~ 6.4 x 10 6 corresponding to ~ 250 hours of livetime. The maximum total exposure is close to the Earth north pole 
with - 2 x 10 9 cm 2 s at 1 GeV. 

In addition to the standard diffuse cuts, corresponding to the post-launch instrument response functions P6V3, we 
applied the following selections: 

• ©nadir < 80° to select photons coming from the earth. 

• ©lat < 65° to avoid the edge of the field of view of the LAT (8lat being the angle with respect to the LAT 
z-axis) . 

III. RESULTS 
A. Intensity Maps 

To generate intensity maps we transform the incident photon directions into the Earth-centric coordinate system 
with angles na dir and Anadir denoting the zenithal and azimuthal angle about the direction pointing from the 
spacecraft towards the center of the Earth, where ^nadir = corresponds to the north. The photons were binned in 
solid angle for energy intervals with spacing 10 bins per decade from 200 MeV up to 500 GeV. The solid angle bins 
were spaced equally in 5° bins in ^nadir with an energy dependent binning in O na dir changing from 1° at 200 MeV 
to 0.1° above 10 GeV to take into account the LAT's energy-dependent angular resolution. The LAT pointing and 
livetime history were evaluated in the same coordinate system in 30-second time steps, creating maps of the exposure 
of the spacecraft to allow conversion from photon counts per solid angle to intensities. 

Figure [Tj shows the resulting exposure and photon intensity maps for 3 different energy bands (left: 200 MeV to 
1 GeV, center: 1 GeV to 10 GeV, right: 10 GeV to 500 GeV). The center of each plot (i.e. 9 nadir = 0°) is the nadir 
direction, while the edge of the plots are the horizon at LAT's altitude (O na dir = 90°). North is to the top of the 
maps and west to the left. All intensity maps have been background corrected by subtracting the average intensity 
from a ring between ©nadir = 79° to 80° averaged over $ na dir- 

The exposure is a factor ~ 2 higher in the north compared to south (top row). This is due to the observation 
strategy in which the majority of the data were taken during a north-orbital pole pointing. The exposure increases 
with energy due to the increase in the LAT effective area. It can also be seen that towards the nadir the exposure is 
much lower than towards the rim. To avoid the introduction of large errors from low exposure bins in the spectral 
analysis, areas where the exposure is less than 0.5% of the maximum exposure in the corresponding energy band are 
excluded from the analysis in all energy bands. The discarded areas due to low exposure match the black regions in 
the center of the exposure maps. Compared with EGRET, the LAT exposure is about a factor of ~ 30 larger fill ]. 

The 7-ray intensity maps (bottom row) in units of photons cm~ 2 s^ 1 sr _1 show the bright limb of the Earth, with 
the emission becoming less intense toward nadir. The Earth limb is very bright with intensities of up to 0.01 photons 
cm -2 s _1 sr" 1 in the low energy bands - more than one order of magnitude brighter than the center of the Galaxy 
for the same energy range. The narrowing of the emission with increasing energy is an effect of both the PSF of the 
LAT improving at higher energies [l5| and the increasing contribution by 7-rays produced in the forward direction of 
the cosmic-ray showering at higher energies. Only 27 photons above 10 GeV are detected with angles Onadir < 60°. 
The maximum of the 7-ray emission is detected from a narrow region at © na dir ~ 68°. The intensity maps also show 
significant modulation between east and west with the effect decreasing as the 7-ray energy increases (this is the 
so-called "east- west effect"). Photons arriving from directions excluded from the analysis due to low exposure are not 
shown on the maps. 

Our measurements are in a good qualitative agreement with expectations from standard cosmic-ray air-shower 
phenomenology. 

The east-west effect that can be seen in the lower energy intensity maps is due to the deflection of cosmic rays in the 
Earth's magnetic field. If there was no magnetic field the cosmic ray intensities at the top of the atmosphere would be 
isotropic, and no east-west asymmetry would exist. However, the presence of the magnetic field biases the cosmic-ray 
intensities interacting in the atmosphere. There are fewer cosmic rays interacting from an easterly direction than the 
west because the magnetic field disallows trajectories from the east. This effect is energy dependent so that low-energy 
cosmic rays are affected more than high-energy cosmic rays (those with energies above ~ 50 GeV are unaffected). 
Since the majority of the emission near the limb comes from viewing the development of the cosmic-ray showers in the 
forward direction, the reduced cosmic-ray intensity from the cast effectively suppresses the shower products coming 
from this direction. This effect vanishes when the atmospheric photons are predominantly generated by cosmic rays 
that have large enough energies to allow them to reach the atmosphere undeflected by Earth's magnetic field. This is 
the case for secondary photons with energies > 10 GeV. 
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FIG. 1: Two-dimensional maps of exposure (top row) and intensity (bottom row) for three different energy bands. The 
coordinate system in these maps is a polar representation such that the nadir-direction (i.e. looking down towards the center 
of the Earth) with lla dir = 0° is at the center of the plot. Onadir is increasing in radial direction with the edge of the plot 
towards the horizon at Gnadir = 90°. The north pole is to the top and west is to the left. The exposure maps are given in units 
of cm 2 s, the intensity maps in units of cm -2 s" 1 sr -1 . 



The high-energy 7-rays {E 1 > 10 GeV) are concentrated in a very narrow band toward, but not exactly aligned with, 
the Earth-limb angle, which is at 66.7° for the LAT's orbital altitude of 565 km. At these high energies the number 
of back-scattered 7-rays formed in the shower development is negligible and the intensity is dominated by forward- 
emitted 7-rays from incident cosmic rays with tangential directions to the Earth's surface. Vertically incident showers 
have their first interaction in the range between 10-20 km, which corresponds to na dir ~ 67°. The atmospheric 
column density in the limb direction is large, being many attenuation lengths of material, which suppresses the 7-ray 
emission close to the limb. In addition, the LAT views mostly grazing-incidence cosmic rays in the higher parts of 
the atmosphere. Therefore the maximum of the emission is in the range between na dir = 67° to 69°. For zenith 
angles larger than about Onadir = 70° the atmospheric depth left to produce a hadronic shower becomes too low 
and, therefore, very few 7-rays are detected. For zenith-angles smaller than the limb angle at na dir ~ 66.7°, the 
Earth and its atmosphere form an opaque barrier to forward-emitted 7-rays. Only the significantly lower-intensity 
back-scattered and off-axis emission from the cosmic-ray shower development is detected for these viewing angles. 
Thus, the maps above 10 GeV show a very narrow band close to but not at the Earth's limb. 



B. Spectrum 

Figure [5] shows the LAT- measured energy spectrum of the atmospheric 7-rays for different nadir bands integrated 
over azimuth angles. The spectra have been corrected for spill-over effects due to the LAT PSF. This was done by 
measuring the intensity in each energy bin outside of the Earth limb, where very little emission is expected due to 
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FIG. 2: Differential energy spectrum for the 7-rays produced in the Earth's atmosphere taken in different O na di r -bands. The 
spectrum of the limb (open circles) shows a power-law behavior. Also shown are 7-ray spectra for the inner parts of the 
Earth's disk (filled circles and points). Indicated are both statistical and systematic errors (grey band), the latter derived from 
studies of the Vela pulsar. The inset shows a comparison to previous data by the SAS-2 satellite [TJ. To allow for a direction 
comparison, the limb-spectrum is also shown for a larger integration region (58° < Snadir < 78°) corresponding to the SAS-2 
measurement. The bottom panel shows the limb-spectrum multiplied by energy to the power 2.75. 



the significantly lower atmospheric density. Outside of the rim is defined here as 70° < O na dir < 73°. Only for the 
low energy part of the spectrum (below ~ 1.5 GeV) this spill-over effect is relevant. Beyond that energy the intensity 
outside the rim is lower by many orders of magnitude than the intensity from inside the rim. 

The average differential spectrum of the limb region (open circles), i.e. for 66° < ©nadir < 70°, shows evidence 
for a turn-over at energies below 3 GeV. Between 3 GeV and 500 GeV the spectrum can be fitted with a power-law 
dN/dE = I E- T with normalization I = (7.4 ± 1.0) x 10 2 MeV" 1 cm" 2 s" 1 sr" 1 and photon index T = 2.79 ±0.06. 
Since both the intensity and the exposure vary considerably over the limb region, a derivation of the statistical 
errors is not straightforward. To determine the statistical errors we use a Monte Carlo technique: 1000 pseudo-count 
experiments are generated for each energy bin, using the LAT exposure and a model of the gamma-ray intensity 
emitted by the Earth albedo derived from EGRET observations [14| . In this model, the average intensity from the 
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limb is renormalized to match the intensities reported here. The RMS of the intensities derived from these pseudo- 
experiments is then used as the uncertainty of the reported gamma-ray intensity from the Earth limb. 

For the fit, the systematic errors (shown as grey band) have been taken into account, conservatively adding them 
linearly to the statistical errors. These systematic errors have been derived from Vela observations by comparing 
the measured 7-rays to various simulation of 7-rays from Vela where these simulations take the uncertainties in the 
instrument response functions into account [la ]. These systematic errors therefore essentially show the contribution 
from the uncertainties in the effective area of the instrument. For energies above 10 GeV, the systematic errors are 
less well known (due to the Vela cutoff) and have been extrapolated from the value at 10 GeV. For this high-statistics 
data set, the systematic errors dominate over the statistical errors nearly up to the highest energies. 

The botto m p anel of Fig. [5] shows the differential intensity multiplied by E 2 - 75 - the spectral index of cosmic-ray 
protons [13, tUf. The x 2 of the fit is 1.2 for 19 degrees of freedom when using the systematic errors in the fitting 
(x 2 /ndf = 89/19 when using only statistical errors). 

The integrated intensity of Earth's limb above 100 MeV is 5.2 x 10~ 3 cm~ 2 s _1 sr^ 1 , which makes the Earth 
by far the single brightest source in the LAT energy range. For this data set there are 16 photons above 500 GeV 
and 3 events above 1 TeV with the highest energy photon having an energy of 1.14 TeV. Events above 500 GeV 
have not been used in the spectral analysis since the energy calibration of the LAT is still under investigation 
at these energies. Also shown as solid circles and dots are the spectra for two inner regions of the Earth's disk 
(60° < 6 na dir < 66° and 30° < O na dir < 60°). The spectra become much softer above a few hundred MeV when 
moving towards the inner part of the disk. Fitting the inner parts by a power-law with energy-dependent index of the 
form dN/dE = I £-r-/3*o S10 (s/M e V) yie id s r = 7J6±3.1 and = 1.78±0.5 for the range between 60° < e nadi r < 66° 
and r = 6.9 ± 0.1 and /3 = 1.90 ± 0.02 for the range between 30° < 6 nadir < 660°. 

The spectrum of the Earth's limb is dominated by 7-rays from cosmic ray interactions in the upper atmospheric 
layers (in tangential directions) pointed towards the LAT that do not suffer large energy losses in the atmosphere (as 
compared to the back-scattered 7-rays). At high 7-ray energies (above 3.2 GeV) the spectrum enters the regime of 
cosmic-ray primaries (> 10 GeV) unaffected by the Earth's magnetic fields, and, therefore, should have a spectral 
index close to that of cosmic rays [l9[ . The fitted power-law index of the atmospheric 7-ray emission from the limb is 
r = 2.79 ± 0.05, compared to an index of T = 2.75 ± 0.03 for the primary cosmic-ray spectrum. Since the spectrum of 
cosmic rays is well measured, eventually, the 7-rays from the Earth's atmosphere might help to calibrate the effective 
area of the LAT beyond the energy regime where we have test-beam data. As previously described, moving from the 
limb to the inner regions of the Earth, the LAT moves from measuring the CR showers from a predominantly forward 
direction to a side-on view and then to a backward direction. This results in a much steeper spectrum for the inner 
parts with very few 7-rays above 1 GeV as shown by the filled circles. 

The inset in the top panel of Fig. [5] shows a comparison with previous results as published in Thompson et al. 
[l3l ]. The spectra are qualitatively similar when taking into account the region for which the SAS-2 spectrum had 
been derived (±10° around the Earth limb). The differences in the measured flux could be due to the limited energy 
resolution of SAS-2 or the fact that when the PSF is folded into the integration region, the match in areas over which 
the measurements were taken is not perfect. 

C. Azimuthal and Zenithal profiles 

Figure [3] shows an azimuthal profile of the limb emission (i.e. for 60° < na dir < 75°) for four different energy 
bands (0.2-1 GeV, 1-10 GeV, 10-30 GeV and above 30 GeV). As can be seen, the east- west effect is a very prominent 
modulation at low energies with a factor of ^ 5 intensity difference between the minimum in the east (at $ na dir = 90°) 
and the maximum in the west (at ^nadir = 270°). This factor is consistent with what has been shown for EGRET 
data [3]. The difference in intensity between the east and the west for the highest energy band is consistent within 
errors with no modulation ((8.0 ±2.0) x 10~ 7 cm _2 s~ 1 sr _1 from the east and (10.0 ±2.0) x 10 _7 cm _2 s _1 sr~ 1 from the 
west). The no- modulation hypothesis provides an acceptable fit only for the > 30 GeV band (\ 2 — 74 for 71 degrees 
of freedom). All other bands show the variation in intensity with azimuth and hence a poor fit to a constant (for 
example, the 10 — 30 GeV band has a \ 2 — 188 for 71 degrees of freedom). A comparison between the intensity from 
the north and the intensity from the south is consistent within errors (variation is below 15%). The behavior of these 
azimuthal variations can be qualitatively understood taking into account the production mechanism. Low-energy 
cosmic-ray shower products are altered by Earth's magnetic field and due to their predominantly positive charge 
preferentially selects particles coming from the west. Higher-energy cosmic rays have higher magnetic rigidities and 
are therefore less affected by Earth's magnetic field. 

Fig.0]shows the zenith angle distribution for several energy bands separated into angular segments of 90° referenced 
to the nadir centred on the directions north, east, south, and west. These plots have been background subtracted by 
taking the average intensity in the wedge at 9 na di r = 80° (beyond the Earth's disk) as a background estimation. For 
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FIG. 3: Photon intensity as a function of azimuth for 4 different energy bands along the limb (60° < lla dir < 75°). Systematic 
errors as derived from LAT data on the Vela pulsar are shown as grey bands. A modulation is apparent in the low-energy 
bands which fades away at higher energies (see text for more details). 



comparison the dashed line shows the LAT PSF (averaged between front and back-converting events) calculated for 
the lower bound of the energy bin (the low energy photons dominate the emission in each of the energy bins). Note 
that for the two higher-energy bins the range and the bin size of these histograms have been reduced to feature the 
limb emission. As can be seen, the Earth emission is wider than the LAT PSF even for the highest band, in which 
the LAT PSF is better than 0.1° (68% containment radius). The good agreement in intensity from north and south 
can be seen in this figure which provides a systematic check on the true nature of the east- west effect. 

The angular distribution of photons near the rim can be used to compare the 7-ray emission with a simple atmo- 
spheric model to investigate the effect of the absorption of 7-rays. The proton-air interaction length is Xi ~ 85 g 



cm~ 2 , which can be estimated using the value ~275 mb for the total inelastic p 14 N cross section [20|. In the thin 
target regime, when the atmospheric depth is X <C Xi, the intensity of produced 7-rays is oc X. A competing process 
is the attenuation of 7-rays due to the e + e~-pair production, with the corresponding scale, the radiation length, of 
Xo ~ 38 g cm -2 in the atmosphere [2lj . The interplay between these processes sets the value of column density when 
the effect of absorption of 7-rays becomes noticeable > 10 g cm -2 . 

Indeed, this behavior is demonstrated in Figure [3 which shows the integrated atmospheric column depth along the 
line of sight for a given nadir angle at a height 560 km for two different atmospheric models. The solid line shows the 
curve for a simple model of the atmosphere with an exponential density profile with scale height 6.8 km [22] following 
the barometer formula X(h) = X(h — Q)e~( h ' h "> with X(h) in units of g cm -2 denoting the column density of air 
overlaying a point at altitude h in cm. The dashed line shows a more realistic atmospheric model (NRLMSISE-00) 
as given by Picone et al. [23j |. Overlaid on the atmospheric column densities is the 7-ray intensity (scaled to match 
the barometric atmosphere profile) averaged over all theta angles (solid circles). Also shown is the deconvolved LAT 
intensity (open circles) in these angles for the PSF at the median energy of this data set (4.7 GeV). The plot shows, 
that far enough out beyond the rim (in the plot for angles 9 na dir > 68.3°), the 7-ray intensity follows closely the 
column density. This region is the thin-target regime, in which the attenuation of 7-rays is not important yet and the 
intensity scales directly with the amount of atmosphere for the proton-air interaction. For column densities larger 
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FIG. 4: Average intensity as a function of nadir angle in wedges towards the different directions (north, west, south, east) for 
4 different energy bands. The narrowing of the emission with increasing energy is seen as well as the change in shape (with a 
broad tail towards the center of the Earth for low energies) and a very narrow limb for the high energies. The comparison of 
the PSF calculated for the lower bound of the energy bin (averaged between front and back-converting events) is shown as a 
dashed line. Note that the range of the histogram and the binning is different for the two high-energy bands. 

than ~ 10 g cm -2 the atmosphere starts to be optically thick for 7-rays and absorption (i.e. a deviation from the 
scaling with density) sets in. 

If the target is thin enough that the secondary 7-rays are practically not attenuated, the observed 7-ray intensity 
should have a power-law index close to the index of ambient cosmic rays. Indeed, it has been shown that the 
assumption of a constant fraction n ~ 0.17 of energy of the incident proton released in the secondary 7-rays (5- 
function approximation) works well in proton-proton interactions for 1 GeV < E~ < 100 GeV [H, The Earth's 
atmosphere consists predominantly of N and O, and there is a significant fraction of He in cosmic rays. This difference 
in the beam and target composition could be taken into account by introducing an approximate energy independent 
correction factor, but we are currently interested in the spectral slope and not in the absolute normalization of the 
7-ray intensity. Therefore, a comparison of the ambient spectrum of protons, with the energy scaled by a factor of 
k, with the observed spectrum of 7-rays may be used to check the equality of their spectral indices and to test the 
energy calibration of the LAT instrument. 

Figure [5] shows the 7-ray spectrum for the optically thin region (68.6° < na dir < 69.6°) compared to the proton 
spectrum with the energy scaled by a factor of k. The proton intensity is scaled by first scaling the energy of each 
intensity point by a factor k, and then re-normalizing the resulting proton intensity down by 0.7. In this way the 
proton intensity matches the 7-ray intensity at 1 GeV (i.e., the ratio is normalized to unity at 1 GeV). The general 
good agreement is further demonstrated by the inset in Figure [5] which shows the ratio of the 7-ray intensity to the 
scaled proton intensity. 
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FIG. 5: Comparison between the integrated column density with the LAT-detected gamma-ray intensity as a function of nadir 
angle for events above 3.6 GeV. Two different model atmospheric profiles have been used: a simple barometric atmosphere 
with scale height of 6.8 km [22l ] shown as solid line and a more realistic atmospheric model from [23fl . Shown are the measured 
gamma-ray intensity (solid circles) and the gamma-ray intensity deconvolved by the PSF of the median energy (4.7 GeV) of 
this data set (open circles). A good correspondence between the 7-ray intensity and the column density can be seen at angles 
©Nadir > 68.3°. For these angles the atmosphere is thin enough (column density < 10 g cm -2 ) so that no significant attenuation 
occurs and the 7-rays are directly related to the amount of target material for the incoming cosmic rays. The absolute level of 
the 7-ray intensity has been scaled to match the column density for the barometric atmosphere model. 



IV. CONCLUSION 



Data from the LAT have been used to study the 7-ray emission generated in the interactions of cosmic rays with 
the Earth's atmosphere for over more than four orders of magnitude between 100 MeV and 1 TeV with unprecedented 
precision. The data set contains 218 photons above 100 GeV, and 16 photons above 500 GeV. These measurements 
therefore demonstrate the capability of the LAT to detect and determine energies for photons up to TeV energies. 
Two-dimensional intensity maps as well as azimuthal and zenithal profiles and energy spectra of the Earth albedo 
emission have been derived which show several effects: 

• A bright limb at the Earth's horizon. These limb 7-rays are dominantly generated by grazing-incidence cosmic- 
ray showers coming directly towards the LAT. The 7-ray spectral index follows the spectral index of the incident 
cosmic ray spectrum in the limb region up to 7-ray energies of E 100 GeV. 

• A soft-spectrum nadir region dominated by 7-rays back-scattered at large angles originating from cosmic-ray 
showers developing deep into the atmosphere. Due to increasing collimation of secondary particles with increas- 
ing energies, these 7-rays show a much softer spectrum. 

• An east-west modulation for energies below ~ 10 GeV, caused by the deflection of primary cosmic rays in 
Earth's magnetic field. Higher energy cosmic rays are less affected by the magnetic field, resulting in a fading 
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FIG. 6: A comparison of the measured gamma-ray intensity for the angular interval 68.6° < G na dir < 69.6° (thin target case, 
gamma- ray attenuation is insignificant) with the scaled cosmic ray proton intensity (see text for details). The inset shows the 
ratio between the 7-ray data and the scaled proton data and shows the general good match between the two measurements. 
The data have been taken from Sanuki et al. [18} and Ryan et al. fl7j . 



modulation with increasing energy. The north-south ratio is equal to 1 within errors. 

Figure [7] shows a compilation of photons generated in the atmosphere of the Earth in the energy range between 
10 keV and 1 TeV. A direct comparison of the measurements from different instruments is not straightforward, since 
different instruments are integrating over different regions of the Earth. Also, the measurements have been taken at 
different times during the solar cycle and therefore correspond to different levels of solar modulation of the primary 
cosmic rays. The peak of the Earth 7-ray albedo energy flux is in the LAT range (for the limb, but potentially also 
when averaging over the whole Earth). The position of this peak can be understood from the power-law index of the 
cosmic-ray spectrum and the fact that neutral pion production is the dominant process. It is similarly visible in the 
spectrum of the diffuse Galactic emission [30J . Evident in Figure [7] is the huge energy range over which the LAT can 
make spectral measurements. 

The LAT data provide a picture of the 7-ray emission from the Earth's atmosphere unprecedented in energy 
range, resolution and statistical precision. The data can be used to understand and model the interaction of 
cosmic rays in Earth's atmosphere and magnetic field. Since the Earth is such a bright source for the LAT, these 
data also provide valuable information for 7-ray background studies and can eventually be used as a calibration source. 
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